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Abstract 

The spectroscopic properties of an open large Bunimovich cavity are studied numerically in the 
framework of the effective Hamiltonian formalism. The cavity is opened by attaching two leads to 
it in four different ways. In some cases, the transmission takes place via standing waves with an 
intensity that closely follows the profile of the resonances. In other cases, short-lived and long-lived 
resonance states coexist. The short-lived states cause traveling waves in the transmission while 
the long-lived ones generate superposed fluctuations. The traveling waves oscillate as a function of 
energy. They are not localized in the interior of the large chaotic cavity. In all considered cases, 
the phase rigidity fluctuates with energy. It is mostly near to its maximum value and agrees well 
with the theoretical value for the two-channel case. 
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I. INTRODUCTION 



During last years, much interest is devoted to the study of the spectral properties of 
small cavities coupled to a small number of channels. As a function of the coupling strength 
between cavity and attached leads, the results obtained show the resonance trapping effect 
in theoretical ^ as well as in experimental [2 1 studies. Short-lived whispering gallery modes 
are formed by attaching the leads in a suitable manner to cavities of different shape with 
convex boundary . These direct processes can result in deterministic transport as signified 
by a striking system-specific suppression of shot noise jj]. The corresponding pathways are 
localized inside the cavities: the whispering gallery modes near to the convex boundary and 
the bouncing ball modes around the shortest pathway between the two attached leads 0. 

The whispering gallery modes can be characterized well by classical values such as path- 
way length and travelling time j^. The last value is related to the lifetime of the corre- 
sponding state, for example to the lifetime of the whispering gallery mode. That means 
the shorter the pathways of the direct modes, the shorter the lifetimes of the corresponding 
resonance states and the better these states are separated from the other resonance states of 
the system by their lifetimes. It follows immediately from these relations that it is easy to 
identify the whispering gallery modes inside a small system, i.e. in a cavity whose attached 
leads have a width that is relatively large as compared to the area of the cavity. It might 
however be difficult to identify them in a large system that is relatively weakly coupled to 
the leads due to its large area and relatively small width of the attached leads. Here the 
pathway between the two attached leads is long. The question arises therefore whether or 
not whispering gallery modes can be identified also in large cavities. 

Guided by the results of previous calculations for small cavities 0, Q], we have chosen 
billiards of Bunimovich type for the study of this question. We opened these billiards by 
attaching two leads to them and used four different geometries for the positions and orienta- 
tions of the leads with the aim to study the transmission through the cavity under different 
conditions. The spectroscopic study is based on the method of the effective non-hermitian 
Hamilton operator i^efr that describes the spectroscopic properties of an open quantum sys- 
tem, i.e. of a quantum system that is opened by embedding it into a common continuum 
of scattering wave functions P]. This means in the present case that leads are attached 
to the closed cavity for the propagation of the scattering wave functions corresponding to 
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one channel in each lead. In some of the open cavities, the expectation value {t{E)) of the 
transmission amplitude t{E) shows an oscillatory behaviour. The oscillation length depends 
on the geometry of the attached leads. It may be large for bouncing ball modes when the 
distance between the input and output leads is small. It is however much smaller in the 
case of the whispering gallery modes that appear in a small attached half stadium and are 
characterized by a relatively large distance between the input and output leads. 

The paper is organized in the following manner. In Sec. II, we provide a few of the 
basic equations that describe the different time scales in open quantum systems. The time 
scales are determined by the lifetimes of the resonance states which are obtained from 
the eigenvalues of the effective Hamiltonian H^s that describes the open quantum system. 
They depend on the manner the leads are attached to the quantum billiard. We show 
further the results of numerical simulations for the transmission through the four different 
cavities. In the case of whispering gallery modes and bouncing ball modes between the 
two attached leads, {t{E)) oscillates with a period that is determined by the momentum 
k and the geometry of the open cavity. In Sec. Ill, we consider the eigenfunctions of the 
effective Hamiltonian. The eigenfunctions are biorthogonal with the consequence that, in 
the regime of overlapping resonances, the real and imaginary parts of the eigenfunctions 
may decouple, to a certain degree. We relate this decoupling to the phase rigidity of the 
scattering wave function which expresses the distortion of the scattering wave function by 
the overlapping of the different resonance states. As for the transmission, we calculate the 
expectation value of the phase rigidity. For isolated resonances, the real and imaginary parts 
of the eigenfunctions are related to one another in the standard manner and the transmission 
takes place by standing waves. In the overlapping regime, however, the real and imaginary 
parts of the eigenfunctions decouple from one another and eventually traveling waves arise. 
When fully evolved, the traveling waves are superposed by long-lived standing waves. The 
traveling waves are described by the so-called optical S matrix. The results are summarized 
in the last section. 
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II. EIGENVALUES OF THE EFFECTIVE HAMILTONIAN AND SEPARATED 
TIME SCALES 



A. Eigenvalues and transmission 

The energies and widths of the resonance states of an open quantum system can be ob- 
tained from the poles of the S matrix or directly from the eigenvalues Zk of the corresponding 
effective Hamilton operator i^gfr fl- For a quantum billiard with two attached leads, the 
effective Hamiltonian is 



H,^ = Hb+ E j,^ ^ rr VCB (1) 

C=R,L ^ - 



where Hb is the Hamiltonian of the closed quantum billiard, He is the Hamiltonian of the 
left (C = L) and right (C = R) lead and = E + iO. The second term of H^s takes into 
account the coupling of the eigenstates of Hb via the modes propagating in the leads when 
the system is opened. It introduces correlations between the states of an open quantum 
system which appear additionally to those of the closed system j^. The (real) eigenvalues 
of the Hamiltonian Hb are the energies of the discrete states of the closed system, 
while the (complex) eigenvalues zx of H^s provide the positions Ex and widths Tx of the 
corresponding resonance states of the open system. There is a one-to-one correspondence 
between the number of eigenstates of Hb and that of Hcs- 

Since the effective Hamiltonian depends explicitely on the energy E, so do its eigen- 
values Zx- The energy dependence is small, as a rule, in an energy interval that is determined 
by the width of the resonance state. The solutions of the fixed-point equations 

Ex = Re{zx)iE=E, (2) 

and of 

Tx = -2lm{zx)\E=E:, (3) 

are numbers that coincide approximately with the poles of the 5" matrix. The width Tx 
determines the time scale characteristic of the resonance state A. The amplitude for the 
transmission in the one-channel case is 0| 

A ^ - -A 



where the eigenf unctions of Hcs are denoted by (px and the scattering wave functions in 
the leads by According to the transmission is resonant in relation to the effective 
Hamiltonian Hcs- This holds for narrow resonance states as well as for the short-lived and 
long-lived resonance states that appear after redistribution of the spectroscopic properties 
of the system. Such a redistribution is studied first numerically in a nucleus 0| and then 
analytically, by using statistical assumptions and neglecting the energy dependence of the 
z., in a large chaotic system 0. It is caused by branch points in the complex energy plane 
9|. That means, the eigenvalues zx of i^efr determine the time scale of the transmission. 

The coupling matrix elements Vbc, Yh b b etween billiard and attached leads can be 
calculated in the tight-binding approach p, |lO|. When they are small, it is Ex ~ and 



A ~ where E^ is the position of the isolated resonance state and r° its width that is 
determined by the Vbc, Vcb- For large Vbc, Vcb, however, E^ and Ex as well as r° and 
Tx may be very different from one another due to reordering processes taking place in the 
system at strong coupling to the environment (full opening of the quantum billiard). In 
this regime, short-lived and long-lived resonance states coexist P] . Examples are short-lived 
bouncing ball modes or whispering gallery modes that may coexist with long-lived resonance 
states in a small quantum billiard 3|. 

Whispering gallery modes may appear in small chaotic as well as regular billiards with 
convex boundary when fully opened and the leads are attached to them in a suitable man- 
ner. Examples are billiards of Bunimovich and circular type j^. The whispering gallery 
modes are localized near to the boundary of the billiard. They have an approximately equal 
distance in momentum k from one another, and their positions in energy are determined by 
the number of nodes which increases with increasing energy. The widths are proportional 
to the length of the pathway along the convex boundary (except for threshold effects) ^. A 
shot-noise analysis has shown that they support direct transport processes j^. They deter- 
mine, therefore, the optical S matrix. The long-lived states however feature indeterministic 
processes corresponding to the universal prediction of random matrix theory [J]. They cause 
the fluctuations of the transmission probability. 
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B. Numerical simulation 



In this section, we show the results of numerical simulations for the transmission through 
a cavity of Bunimovich type with leads attached in four different ways. The calculations for 
the transmission are performed by using the boundary element method ^|. The eigenvalues 
of Hpff are obtained by applying the tight-binding lattice Green function method given in 



Ref. 



12l | and using the general relation between H^s and the Green function. The energies 



and widths of the resonance states are obtained by solving the fixed-point equations (j21) 
and 0. The ensemble average (t) is performed from 200 different positions of an internal 
obstacle by keeping fixed the area of the billiard. The energy interval considered is divided 
usually into 20 energy bins. 

In Fig. we show the ensemble expectation values {t{E)) of the transmission as a 
function of the energy E in the region E = 14 — 26 for the four cavities shown in the 
insets of Fig. |21 The oscillating contribution from the short-lived states as well as the shifts 
between Iie{t{E)) and lm{t{E)) can be seen clearly in the cavities (c) and (d). Interesting is 
the geometry of the cavity (b) where the transmission {t{E)) changes its nature at E ^ 21. 

In Fig. 121 we show the eigenvalues of the effective Hamiltonian for the four cavities. 
The values Ex and Tx of the short-lived states are calculated by solving the fixed-point 
equations (j21) and 0. We have resonance states with well separated time scales in Figs. 
|21^c) and (d). In these cases, {t{E)) is large and oscillates. In Fig. Efa), separated time scales 
can not be identified and {t{E)) is relatively small. In Fig. Efb), separated time scales can 
be identified but the difference between the short-lived and the long-lived states is smaller 
than in (c) and (d). Furthermore, the widths of the long-lived states are spread in (b) over a 
comparably large range and the widths of the short-lived states show an irregularity around 
E ^ 21. At this energy, {t{E)) changes its nature as can be seen from Fig. ^ (b). In any 
case. Figs. H^b) and|2fb) show that the sensitivity of {t{E)) against parameter variations 
is large when the direct pathway between the input and output leads is large and not well 
separated from other pathways through the interior of the cavity. This sensitivity can be 
seen also in comparing the results for the cavities (a) and (d). While there is almost no 
separation of the pathways through the small attached half stadium from those through the 
large Bunimovich stadium in (a), both parts are well separated in (d). As a consequence, 
we see whispering gallery modes in the small attached stadium in (d) but not in (a). 
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FIG. 1: The ensemble expectation value {t{E)) as a function of the energy E for the cavities (a), 
(b), (c) and (d) shown in the insets of Fig. [2 The cavities consist of a Bunimovich stadium 
connected to two waveguides directly, as in panels (a), (b) and (c), or through a smaller half 
stadium, as in (d). Full lines: Jie{t{E)), dashed lines: lm{t{E)). Most {t{E)) show oscillations 
that are related to the distances between the input and output leads. The energy is in units /2m 
and d = 1 is the width of the wave guide. 

The short-lived states determine the value of the so-called optical S matrix. 
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FIG. 2: The solutions of the fixed-point equations @ and (j31) for the resonance states of the four 
different cavities. The short-hved states are marked by crosses, the long-hved ones by dots. A clear 
separation of the time scales can be seen in (c) and (d). The neighboring short-lived resonance 
state in (c) lies at Ex — i/2 Tx = 27.74 — 6.89 i. In the insets, the cavities are shown. In order 
to see the differences between the four open cavities, the attached leads are also shown up to an 
arbitrary distance L. The eigenvalues are calculated with L = 0. 

III. EIGENFUNCTIONS OF THE EFFECTIVE HAMILTONIAN AND PHASE 
RIGIDITY 

A. Eigenfunctions and transmission 



The scattering wave function is solution of the SchrodinErer equation (H — E)"^ 
in the total function space with the hermitian Hamilton operator H. It reads [5|, |6|] 



c 



E 

A 



E+ -H, 



c 



E-zx 



(5) 
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The (px are complex and biorthogonal 



(0a|0A') = (0II0A') = ^A,A' (6) 

\{<t>x\<l>x)\ = A,>l; \{(t>x\(t>x')\ = B^' > . (7) 

Eq. shows that the scattering wave function in the interior region of the quantum 
dot is determined, above all, by the complex eigenf unctions 4>\ of i/efr, 

*^(r)-E0AW^^^- (8) 

A ^-^^ 

At the energy E ^ E\, the eigenfunction (f)\{r) of the effective Hamiltonian Hcs gives the 
main contribution. This holds true especially at the energy of a narrow resonance state. As 
a numerical example, the eigenfunction (px of a whispering gallery mode is shown in Fig. El 
for the cavity (d) and compared with the corresponding scattering wave function at the 
energy Ex- The scattering wave function contains contributions from other eigenf unctions 
also at E = Ex- Nevertheless, it shows, at this energy, the typical structure of the whispering 
gallery mode. As can be seen from the figure, it is localized not inside the large Bunimovich 
cavity but outside of it in the small half stadium. 

The real and imaginary parts of the eigenfunctions <px are more or less decoupled in the 
regime of overlapping resonances The value 



^A 



Jdr{\ReMr)\'-\^^Mr)\y 



^x\9x) 



(0a|0a) 

mm 



Jdr{\ReMrW + \^^MrW 
is a measure for the biorthogonality of the eigenfunctions of Hpf^. The phase rigidity \p\'^ of 
the scattering wave function \l/(r) is considered in e.g. Refs. 



^ /cir vl>(r)^ ^ (|Re§(r)p-|M(r)p) 
^ /rfr|^(r)|2 ^ /(ir (|Re^^(r)|2 + |M(r)|2) ^ ' 

where 6 is an angle providing that Re\^(r) and Im\E'(r) are orthogonal. The value p is related 
to the rx according to (jH)). 

For an isolated resonance state, ^4^ ~ 1 and ~ 1 at the energy E = Ex- At this energy, 
the transmission probability shows a peak. Approaching a branch point in the complex 
energy plane where two eigenvalues and zx^ coalesce, — > oo and rx ^ jlsHlot. 
Here, the widths bifurcate : one of the states aligns with the channel wave function and 
becomes short-lived while the other one becomes long-lived jsj. Eventually, the short-lived 
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FIG. 3: The eigenfunction of the effective Hamiltonian i^cff at the energy E = E\ of the 
whispering gallery mode A (left) and the scattering wave function at the same energy E 
(right). E\ — i/2 T\ = 19.53 — 2.25 i. The (j)x is shown up to the attached lead {L = 0) while the 
is shown, for illustration, also in the lead up to an arbitrary finite value L. 

and long-lived resonance states differ strongly from one another and do not cross in the 
complex energy plane . Therefore again, as for non-overlapping resonances, ^ 1 and 
ta — 1 at the energies E = Ex of the long-lived states. In the transmission through the 
cavity, the short-lived states determine the smooth 'background' while the long-lived states 
cause the superposed peaks (fluctuations) in the transmission probability 

This picture can be translated to that of 'standing' and 'traveling' waves. Standing waves 

N 

tP{r) = (2iV)"^/2 j2 cos [{On + kn ■ r)] (11) 



n=l 



cause the Porter-Thomas statistics for the intensity [17|. Tuning the frequency to a reso- 
nance, an intensity pattern is generated that closely follows the profile of this resonance. In 
this situation, the resonances are isolated from one another and ^4^ ^ 1, ^ 1. 

When the cavity is fully open, the local field can be viewed as a sum of a number of 
traveling modes ariving at a point from various scattering processes , 

N 

^{r) = (2iV)-i/2 ^ exp [z(^„ + K ■ r)] (12) 

n=l 

where the phases 6'„ are completely random and the wave vectors A;„ are uniformly dis- 
tributed. Both Re('i/') and lm('0) are independent Gaussian variables what leads to the 
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Rayleigh distribution for the intensity /(r) = [-(/^(r)!^. It applies to a monochromatic wave 
propagating in an open system ('travehng wave' excited by a monochromatic source jl7|). 
In the corresponding description of this situation with the effective Hamiltonian formalism, 
short-lived and long-lived resonance states coexist in the system and determine, respectively, 
the smooth 'background' and the superposed peaks (fluctuations) of the transmission proba- 
bility 0]. The short-lived resonance states are strongly related to the channel wave functions 
(scattering states in the leads) due to the large overlap integral of their wave functions with 
those of the channel wave functions. The transmission induced by these states shows there- 
fore the same dependence on the momentum k as the scattering wave functions in the 
leads ('traveling' waves). Obviously, the 'monochromatic source' by which the traveling wave 
is excited according to is, in the one-channel case, the channel wave function since 
it causes the alignment of one of the wave functions (j)x in approaching the branch point in 
the complex energy plane P]. The traveling waves determine the optical S matrix. 



B. Numerical simulation 
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FIG. 4: The ensemble and energy averaged phase rigidity for resonance states of the cavity 
(c). Energy interval [22, 23]. The full line is calculated from the equation for the phase rigidity 
distribution in the case of a chaotic cavity with energy averaging and 2 channels [ij] . The results 
for the cavities (a), (b) and (d), see insets in Fig. [21 in the same energy interval [22, 23] as well as 
those for the cavity (b) in the energy interval [20.5, 21.5] are nearly the same. 
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The distribution of the phase rigidity is calculated by means of (jlUj) . We take the ex- 
pectation value ( |/9p) for an ensemble of 200 cavities with different positions of the obstacle 
in the interior and for 20 different energy values inside the energy interval considered. The 
results are almost the same in all cases considered, i.e. for the four different open cavities in 
the energy interval [22, 23] and, in addition, for the cavity (b) in the energy interval [20.5, 
21.5]. A typical result is shown in Fig. |3J It a gree s well with the theoretical value (full line 
in Fig. for a chaotic cavity and 2 channels [l5[. The phase rigidity is mostly near to its 
maximum value in the two-channel case. This corresponds to the fact that the transmission 
is caused by either standing waves or traveling waves with superposed fluctuations. The 
transition between the two scenarios takes place in a comparably small region according to 
the examples studied in 0, Q • 



IV. SUMMARY 



The spectral properties of an open cavity depend strongly on the manner the leads are 
attached to it. We studied the eigenvalues and eigenfunctions of the effective Hamiltonian 
Hcs describing a large cavity of Bunimovich type with two leads attached in four different 
ways and one channel in each lead. The transmission is resonant in all cases in relation 
to the effective Hamiltonian of the open quantum system. In some cases, the transmission 
takes place via standing waves in the cavity with an intensity that closely follows the profile 
of the resonances. In other cases, two different types of resonance states appear which differ 
by their lifetimes. The short-lived states cause traveling modes while the long-lived states 
appear as fluctuations of the transmission probability. The short-lived states, including the 
whispering gallery modes, can not be identified in the interior of large cavities where the 
pathway between input and output leads is large. Therefore their widths are relatively small 
in this case, and it is impossible to identify them in the 'sea' of long-lived resonance states. 

The optical S matrix is related to the short-lived states (traveling waves) as can be seen 
from the eigenvalues and eigenfunctions of H^s. This relation is however not necessarily true 
also in the opposite direction : the eigenvalues of i^efr may show separated time scales while 
the optical S matrix is, nevertheless, small. In such a case, the short-lived modes exist, 
according to our numerical results, inside the cavity, and their spectroscopic properties are 
very sensitive against small parameter changes. In the four cases studied by us, ballistic 
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modes do not appear in the interior of the large chaotic cavity. 

In all considered cases, the phase rigidity fluctuates as a function of energy but is mostly 
near to its maximum value. The distribution is characteristic of the two-channel case. 
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